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Abstract 

Direct Sun light absorption by trace gases can be used to quantify 
them and investigate atmospheric chemistry. In such experiments, the 
main optical apparatus is often a grating or a Fourier transform spec- 
trometer. A solar tracker based on motorized rotating mirrors is also 
needed to direct the light along the spectrometer axis, correcting for 
the apparent rotation of the Sun. Calculating the Sun azimuth and 
altitude for a given time and location can be achieved with high ac- 
curacy but different sources of angular offsets appear in practice when 
positioning the mirrors. A feedback on the motors, using a light posi- 
tion sensor closed to the spectrometer is almost always needed. This 
paper aims to gather the main geometrical formulas necessary for the 
use of a widely used kind of solar tracker, based on two 45° mirrors 
in altazimuthal set-up with a light sensor on the spectrometer, and to 
illustrate them with a tracker developed for atmospheric research by 
our group. 

1 Introduction 

Spectroscopic studies of direct Sun light are commonly used in atmospheric 
research. Such experiments make use of the Sun as a light source to quantify 
molecular absorptions in the atmosphere and then retrieve trace gases abun- 
dance. Stratospheric ozone [Ij and greenhouse gases [2] are routinely mea- 
sured with this technique from ground-based Fourier transform infrared (FTIR) 
spectrometers, e.g. within the Network for the Detection of Atmospheric 
Composition Change (NDACC, http://www.ndacc.org/). In the UV-visible 
range, light scattering is more important and enables spectroscopic studies of 
the atmosphere in other geometries such as zenith measurements [3] . Direct 
sun light is however also used 'S] , its unique and unambiguous light path 
making it advantageous for some applications [6j. Beside the spectrometer, 
the main part of the involved apparatus in direct sun light spectroscopy is 
the solar tracker, required to compensate for the Sun's apparent movement. 
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Several kinds of trackers, sometimes referred as heliostats, are used for 
atmospheric spectroscopy, based on set-ups of one or several rotating mir- 
rors. Some of them are equatorially mounted, like in Table Mountain Facility 
[3] or Harestua [7]. In this case, one rotational axis is parallel to the Earth's 
axis. It enables a high tracking accuracy without a computer, since only 
one axis has to be driven at the Earth's rotation speed. To our knowledge, 
it is the only set-up working without retroaction on the Sun's position. On 
the other hand, equatorial mounts are large, need to be aligned accurately 
and their mechanical design is difficult. Most of the trackers used today 
are controlled by a computer which makes possible remote operation and 
automatization. The computer calculates first the Sun position, moves the 
mirrors to point the Sun and then controls this mirror to optimize the signal 
on some kind of light sensor. 




Figure 1: Geometrical set-up of the considered solar tracker, using two 45- 
degree mirrors. Ml and M2, rotating along orthogonal axes. Mirror MO 
directs the Sun light into a spectrometer. A fraction of the light beam is 
deflected toward a 4-quadrant photodiode enabling a closed-loop control of 
the mirrors position. 
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Figure 1 shows a popular altazimuthal tracker design. It consists of two 
elliptical mirrors held in 45 degrees facing each other (Ml and M2). Both 
Ml and M2 rotate along the azimuthal axis and M2 rotates as well along 
the altitude direction. MO, and possibly other fixed mirrors, sets the light 
beam into the spectrometer optical axis. A 4-quadrant photodiode is used 
as a position sensor for a closed-loop control of the mirrors position once the 
Sun's calculated position has been reached with enough accuracy, i.e. once 
the Sun's image is visible by the photodiode. This altazimuthal set-up is 
used with FTIR systems e.g. in Kiruna [8j and Park Falls, Wisconsin [9J, it 
has been installed in Harestua to replace the equatorially mounted system 
|10j . Compact versions have also been developed for field campaigns |lH ll2j. 
A commercial version is sold by Bruker to be installed above their FTIR 
spectrometers [13]. A recent progress has been reported in the pointing 
accuracy [13], replacing the traditional quadrant diode with a CCD camera, 
but the problems discussed thereafter remain the same. 

Because developing a solar tracker is typically a master's thesis work 
[U [T2l lllj . technical implementations are difficult to access in the litera- 
ture. Some more information is available about the systems used in solar 
energy applications but their geometry differ [151 US] • Someone building a 
Sun tracker can quickly find ephemeris calculations in many programming 
languages, but other issues arise quickly. It is first necessary to character- 
ize the field-of-view(FOV) of the 4-quadrant diode in the considered optical 
design. This serves two purposes: determining the accuracy needed for the 
ephemeris's algorithm and making sure this FOV is larger than the Sun's 
apparent diameter (9 mrad). This last point is important to track constantly 
the center of the Sun, which reduces the uncertainties in the air mass factor 
and avoid Doppler shifts on the edges of the Sun ([H]). A second problem 
lies in the correction of the tracker orientation compared to the altazimutal 
system in which the ephemeris is given, necessary for the calculated mode if 
the ground is not leveled. Thirdly, the relation between the quadrant signal 
and the correction to apply on the mirrors position depends on the tracker 
position itself [14j. Understanding this relation is compulsory to achieve a 
smooth tracking. This article deals with these three problems successively. 

2 Theoretical basis 

2.1 Ephemerisis accuracy and field of view of the 4-quadrant 
diode 

Calculating the Sun position in the sky given the time of observation and the 
geographical coordinates is well documented. A reference algorithm with a 
typical error as low as 0.0003° is given by Jean Meeus in [T7], for which 
C(|18]) or Matlab (p^) versions are available. This accuracy degrades with 
higher zenithal angle due to atmospheric refraction, which depends on local 
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meteorological conditions [20]. On the other hand, the absolute accuracy of 
commercial rotation stages used in Sun trackers, e.g Newport RV-160, is only 
0.01°, which reduces the interest of using a complex algorithm. The lowest 
acceptable accuracy is determined by the field of view of the 4-quadrant 
diode: once the Sun image hits the quadrant, the tracking can be performed 
in closed- loop. 



4-quadrant 

mirror M4q |ens L1 diode 




L D f 



Figure 2: Optical scheme from the tracker mirrors to the 4-quadrant diode. 
The field of view seen by the diode depends on the different aperture sizes 
and path lengths. 



Figure 2 shows the typical optical scheme between the tracker and the 
quadrant diode, the optical axis has been aligned for the sake of clarity. Close 
to the spectrometer, a part of the beam from the tracker, whose diameter is 
limited by the mirror diameter $7^ is deflected by a mirror of diameter $ m 
to a lens which focuses the beam onto the quadrant. The maximum field- 
of-view seen by the diode(-FOVi) depends on the focal length of the lens(/) 
and on the diameter of the quadrant (^d), according to FOVi = arctan 
The lens LI is unlikely to reduce the FOV assuming its size superior to 
the mirror's one. Indeed, the beam is parallel before the lens which implies 
that distance D can be reduced if necessary. The tracker aperture is more 
important, especially since the diameter of rotation stage is limited and 
the distance between the tracker and the deflecting mirror(M4) depends 
on the observatory configuration. The mirror's FOV due to the tracker is 
FOV2 = arctan^. 

Considering our set-up in Brussels, which is a typical FTIR station, the 
tracker's mirrors are 10 cm wide and the optical path to the spectrometer 5m 
long, which leads to a FOV2 of 20 mrad. On the other hand, the quadrant 
diameter is 6 mm while the focal length of the lens is 200mm, which gives 30 
mrad for FOVi. The actual FOV is the minimum, 20 mrad, limited by the 
tracker size. This value is superior to the apparent diameter of the Solar disk 
which is important to track the center of the Sun. A simple algorithm can 
achieve such an accuracy for the ephemerisis calculation, like the one given 
by N0AA(|21|). It is however necessary to take into account the orientation 
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of the tracker, which can lead to pointing errors superior to the FOV. 
2.2 Correcting the tracker orientation 

A source of error when pointing to the calculated Sun position originate from 
the orientation of the observatory compared to the altazimuthal system, 
If the ground is not completely leveled, this problem add offsets to the 
calculated positions, which are not constant during the day. To take that 
into account in the calculation requires to determine the Euler Angles of the 
observatory frame compared to the altzimuthal system. We discuss their 
effects before describing our way to determine them. 




South 



Figure 3: Illustration of the Euler Angles of an observatory compared to the 
altazimuthal coordinates system. 

Euler angles represents the consecutive rotation on the three axes needed 
to account for the pitch(7), roll(/3) and yaw(a) of the observatory. Convert- 
ing the solar altitude and azimuth to the observatory frame requires thus 
to compute the multiplication of three rotation matrices along the different 
axes, i.e Rxi'y), Ryil^) and Rz{a)- 
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The calculations leads to: 



^offsets = (2) 

cos a cos P — sin a cos P sin (5 

cosQ;sin;Ssin7 + sinQ;cos(7) cos a cos 7 — sin a sin ^ sin 7 — cos ^ sin 7 

sin a sin 7 — cos a sin P cos 7 cos a sin 7 + sin a sin j3 cos 7 cos ^ cos 7 



In Cartesian coordinates, the position of the SmL{xt,yt,zt) in the obser- 
vatory frame will thus be related to the solar spherical coordinates(a2;o,aZto) 
in the altazimuthal system: 



M, 



offsets 



X 



cos alto COS azQ 
cos alto sin azQ 
sin alto 



(3) 



In the above equation, the distance from the observatory to the Sun, r, 
is set to unity. This does not matter for the application and simplifies the 
calculation. Substituting Moffset with Eq. 2 we get the following expressions 
for those coordinates: 



(4) 



Xf = (cos a cos alto cos azQ — sin a cos alto sin azo) cos P 
+ sin alto sin (3 

yt = (cos a cos alto sin azo + sin a cos a/to cos azg) cos 7 

— (cos Q cos alto cos azo — sin a cos a/io sin azo) sin /3 sin 7 

— sin alto cos /3 sin 7 

= (cos a cos aZto sin azo + sin a cos aZfo cos azo) sin 7 

—(cos a cos aZto cos azo — sin a cos aZto sin azo) sin /3 cos 7 
+ sin aZto cos P cos 7 



These new Cartesian coordinates can then be converted to altitude (a/tt) 
and azimuth(o^;t) angles relative to the tracker: 
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(5) 



In the above equation arctan4Q stands for the variation of the arctangent 
function cahed atan2 in many programming languages. 

Determining Euler angles accurately by measurements is not easy. An 
analytical method to estimate them is given in |15) which basically consists 
in recording the position of the tracker at three different times and solving 
Eq. m This is appropriate for the studied case, i.e. a collector for solar 
energy application installed outside with only one mirror and no closed-loop 
control. With our considered two-mirror tracker, which does not collect 
light but directs it toward a spectrometer, other source of misalignments 
appear. Indeed, the tracker is also likely to be misaligned compared to the 
spectrometer, and mirror themselves can be tilted. Other angles can be 
considered in Moffsets and is done in [8]. However as some of them change 
with time (due to thermal expansion for instance), it is not very useful. In 
practical applications the three offsets are likely to enable the calculated 
mode to reach the Sun within the FOV of the 4-quadrant diode. With the 
closed-loop control it is easy to track the Sun a whole clear-sky day providing 
an operator correctly sets the Sun tracker initially. Euler angles can then 
be fitted all along the recorded positions of the mirrors during the day. It 
has the advantage that other sources of misalignment are considered; even 
if only three angles are fitted which may not exactly be the Euler angles, 
they minimize simultaneously the effects of all offsets. We implement this 
method in Sect. [3J This requires the closed- loop control of the tracker on 
the Sun position. 

2.3 Active tracking 

The photodiode signal indicates how tilted the Sun beam is compared to 
the optical axis of the spectrometer. This information must be converted 
into angular displacements on the altitude and azimuth axes of the tracker 
to correct the misalignment. If the photodiode was placed on the reference 
frame of the mirror M2 this conversion would be straightforward, but due to 
its position after the tracker it depends on the position of the tracker mirrors. 
A try-and-error method to correct the misalignment is theoretically possible 
but a smoother tracking can be achieved if the conversion is understood. 

The conversion can be expressed once again by a matrix, which transform 
in this case a vector hitting mirror Ml to a vector pointing to a direction 
in the sky given by its altitude and azimuth. It is the opposite of the light 
direction but is simpler to figure out and considering Fermat principle, yields 
the same information. 
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Figure 4: The tracker mirrors and their rotation can be modeled as rotation 
matrices in their reference frames, which apply to the beam vector. 



As with the Euler Angles, the transformation reduces conveniently to a 
combination of rotation matrices. For instance, compared to the reference 
frame 5Ro, the reflection on mirror Mq is expressed as i?y(7r/2) while the 
azimuthal rotation of an angle 9i is written Rz{9i)- 

Figure 4 presents the tracker pointing to an azimuth 9i and a zenith 
angle of ^2- The optical system inside the frame, with only mirrors Mi and 
M2 can be expressed as a transformation whose matrix Mtracker is: 

-^tracker 



Rz{ei) X Ry{e2) X Rx{-) x Ry{-92) x Rx{--) x Rz{-9i) (6) 



i.e: 
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Developing the matrix product yields to the matrix of the tracker optical 
system as a function of the tracker position {01,02)'- 

^tracker — 

cos dl cos 92 cos(Si+92)+sin 9\ sin(Si+92) cos 9\ cos 02 sin(9i+92)— sin 9i cos (Si +62) cos 9\ sin 62 

sin 9\ cos 92 cos (Si +62)— cos d\ sin(9i+S2) cos 9\ cos(Si+S2)+sin 9\ cos 62 sin(Si+92) sin 9\ sin 62 
— sin 02 cos 61 +02 — sin 02 sin 01+02 cos 02 

The transformation expressed by Mtracker can now be applied to a vector 
corresponding to the Sun light beam direction on the spectrometer side of 
the tracker, it will lead to the position of the Sun in Cartesian coordinates. 
The vector is built from the diode 4 signals (VA,VB,VC,VD), as represented 
in Fig. 5. Basically an offset position (ei,e2) is computed for the Sun spot 
on the diode plane compared to its center by: 




Figure 5: Sun spot hitting the quadrant, not to scale. 
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ei = [Vb + Vc) - {Va + Vd) 
£2 = {Va + Vb) - {Vc + Vd) 



(8) 



The spot offset (£1,62) defines 2 coordinates of tfie beam vector. Tlie 
last one, A, sfiould represent tlie distance from tlie diode to mirror Ml. 
Multiplying Mtracker by the quadrant vector (ei,e2,A) would yield accurate 
Sun angles after conversion to spherical coordinates, but is not practically 
possible with a diode, contrary to an imaging sensor. The calculated position 
{xs,ys,Zs) hereafter is thus not absolute but is sufficient to get the sign of 
the rotations to apply on the axes. A can be chosen arbitrarily as long 
as it its absolute value is large enough compared to ei and £2- The solar 
pseudo-coordinates are then: 



^tracker ^ 



(9) 



In practice, the quadrant vector may differ from (£1,62, A) due to reflec- 
tions such as on the mirrors MO and M4q on Fig. 1, necessary to deviate a 
part of the beam to the 4-quadrant photodiode. Defining the position vector 
requires thus to pay attention on the optic path from Ml to the photodi- 
ode, but this is rather easy and reduce to a rotation once the image of the 
photodiode under Ml is figured out. In the next section, we explain how we 
deal with the problem in our particular case. 

Developing Eq. [9] yields to: 



Xs = (cos 6*1 cos 62 sin (^i + 62) — sin 61 cos {9i + 6'2))e2 
+(sin 6*1 sin {0i + 6*2) + cos 9i cos 62 cos (^i + 92))£i 
+A cos 9i sin 62 

Us = (sin 6*1 cos6sin(6'i + 62) + cos 6*1 cos (6*1 + 6'2))e2 ^^^-^ 
+{smdi cos 62 cos (^1 + 62) — cos 9i sin (^i + 62))£i 
+A sin 61 sin 62 

Zs = - sin 62 sin (^i + 6'2)e2 - sin 6*2 cos (^i + 62)61 
+A cos 02 

V 

It is then possible to calculate roughly an altitude(02S') and azimuth(0i5) 
for the Sun applying the Cartesian to spherical coordinates conversion (Eq. 
5). This position is approximate and relative to the tracker since it does not 
take into account the Euler angles described in the last section, but what 
matters is the signs of the differences between these calculated values and 
the current altitude and azimuth relative to the tracker, defined by 9i and 
02- The angular corrections to apply on the two axes are then: 
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del = sgn{eis - 6i)ki 
de2 = sgn{92s - 6*2)^2 



(11) 



Where ki and /c2 are the tracking angle steps that should be small to 
have a smooth tracking, yet large enough for the mechanical resolution of 
the rotation stages and the apparent movement of the Sun. The azimuth 
changes for instance at a rate of 15° per hour, assuming 1 second between 
the steps, ki should not be under 0.004°. 

3 Application for a FTIR measurement station 

Our group has been doing FTIR measurements at Reunion Island for sev- 
eral years (|22).|23j,|24j). The place is interesting since atmospheric measure- 
ments are sparse in the tropical and subtropical region. Aiming at long-term 
monitoring and cost-effectiveness, a station at Saint-Denis (20.9°S, 55.5°E, 
50m a.s.l.) has been automatized [25], which includes sun tracking, mete- 
orological logging and FTIR measurements with a Bruker 120M. The solar 
tracker currently used is a commercial one from Denver University. Since 
September 2009, this station is officially part of the NDACC network. We 
plan to install a second FTIR on the Maido Observatory(21.1°S, 55.4°E, 
2200 m a.s.l) which will start operations in early 2012. This station, which 
will also be automatized, will be based on a Bruker 125 HR spectrometer, 
more appropriate to measure CO2 atmospheric loading, in the framework 
of the new Total Carbon Column Observing Network (TCCON). The spec- 
trometer is now tested at Brussels with an altazimuthal Sun tracker built at 
our institute, which was used to validate the methods described in the last 
section. 

The new solar tracker uses a Newport RV-160 rotation stage for the az- 
imuth rotation and a Vextra stepping motor with a gear box for the altitude. 
Both rotations are driven by a Newport XPS controller, linked to the con- 
trolling PC. The tracker mirrors are elliptical with a 10 cm minor axis. The 
photodiode set-up was purchased from Bruker with the spectrometer and 
is installed at the input window of the spectrometer. It consists of a 1 cm 
mirror which reflects a small portion of the incoming light to a 18 cm focal 
length lens which focuses the beam onto the 4-quadrant photodiode. The 
optical path is shown in Fig. 6. The FOV of the 4-quadrant photodiode is 
20 mrad (see Sect. 12. ip . The algorithm used to compute the ephemerisis is 
the one given by NOAA ([21]). During operation, the mirrors positions are 
refreshed every second according to the calculated position or to the signal 
on the 4-quadrant photodiode using the methods described in 12.31 From 
Fig. 6, it is clear than the beam from the tracker undergo two rotations of 
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4-quadrant diode 
spectrometer window 



Figure 6: Optical paths of our particular set-up from the tracker to the 
spectrometer and the 4-quadrant diode. 



90° along the same axis before hitting the photodiode, which we take into 
account inverting Eq. [TTl 

Figure 7 shows a fit of the Euler Angles. The track was performed 
on 8 February 2011 at BIRA-IASB in Brussels. The left panel shows the 
calculated mirror positions neglecting Euler Angles, together with the ac- 
tual one when the active tracker was operational. Weather was clear-sky 
and enabled, despite the short winter daylight, to record a long Sun path, 
demonstrating the capacity of the active tracking algorithm described in 
Sect 12.31 Around one hundred points were extracted from the log file of the 
tracker position and used to fit the Euler Angles with a unconstrained non- 
linear minimization. The pitch, roll and yaw were respectively estimated 
to be -0.128° ,-0.736° and -1.22°. The right panel shows how these fitted 
angles improve the calculated Sun position. The maximum offset between 
calculated and actual position is now 0.5°, i.e 9 mrad, which is under the 
20 mrad of the photodiode's FOV. Providing this accuracy in the calculated 
mode, the tracking system is able to set the sun's image onto the 4-quadrant 
photodiode and then start the active tracking without the need of an oper- 
ator. 
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Figure 7: Fit of the Euler Angles. 



4 Conclusion 

We have derived the geometry formulas needed to track the Sun with a kind 
of altazimuthal tracker widely used in atmospheric remote sensing. The set- 
up is based on two rotating 45° mirrors facing each other and a 4-quadrant 
photodiode involved in a closed-loop control of the tracker. After discussing 
the required accuracy for the calculated mode and calculating the FOV of 
the sensor, we described how to take into account and estimate the Euler 
angles, representing the orientation of the tracker compared to the ground. 
These sections can actually be applied to other tracking set-ups. On the 
other hand, even if the method is general, the formula for the active tracking 
depend strongly on the optical configuration and may not be used for other 
tracker's geometry. Finally, we have tested the formulas with a custom- 
built solar tracker that will be installed together with a FTIR spectrometer 
at Reunion Island in 2011. 

The next steps involve the determination of the actual tracking accuracy 
with, for instance the Doppler shift of the Fraunhofer lines of the Sun. In 
case it is needed, the tracking smoothness could be improved with a PID 
control instead of the basic retroaction scheme presented here. Another im- 
provement, which can not really worked out at our mid-latitude observatory 
in Brussels, would consist in finding a solution for the singularity occurring 
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at zenith. Indeed, when the solar zenith angle becomes to small, the azimuth 
has very little effect and the closed-loop tracking gets unstable. This is a 
common problem with altazimuthal system and leads to a dead measuring 
time of around one hour at local noon for the NDACC station already ex- 
isting at Saint-Denis. We will deal with this problem if it is significant for 
our new tracker. 
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